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We propose a thermodynamically motivated measure of gravitational entropy based 
on the Bel-Robinson tensor, which has a natural interpretation as the effective super- 
£NJ . energy-momentum tensor of free gravitational fields. The specific form of this measure 

differs depending on whether the gravitational field is Coulomb-like or wave-like, 
Q ■ and reduces to the Bekenstein-Hawking value when integrated over the interior of 

a* a Schwarzschild black hole. For scalar perturbations of a Robertson- Walker geometry 

5— i ' we find that the entropy goes like the Hubble weighted anisotropy of the gravitational 

field, and therefore increases as structure formation occurs. This is in keeping with 
our expectations for the behaviour of gravitational entropy in cosmology, and provides 
■ a thermodynamically motivated arrow of time for cosmological solutions of Einstein's 

. field equations. It is also in keeping with Penrose's Weyl curvature hypothesis. 

in 

^ ■ 1. Introduction 

m 

A key question in cosmology is how to define the entropy in gravitational fields. A 
suitable definition already exists for the important case of stationary black holes pQ, 
^ \ but in the cosmological setting a well-motivated and universally agreeable analogue has 
yet to be found. Addressing this deficit is an important problem, as in the presence of 
gravitational interactions the usual statements about matter becoming more and more 
uniform are incorrect. Instead, structure develops spontaneously when gravitational 
attraction dominates the dynamics [2j |3]. This behaviour is crucial to the existence of 
complex structures, and indeed life, in the Universe. The question then arises, how can 
evolution under the gravitational interaction be made compatible with the second law of 
thermodynamics? If the second law is valid in the presence of gravity, such that entropy 
increases monotonically into the future, then the current state of the universe must be 
considered more probable than the initial state, even though it is more structured. For 
this to be true, the gravitational field itself must be carrying entropy. 

For a candidate definition of gravitational entropy to be compatible with 
cosmological processes, such as structure formation in the Universe, it needs to be 
valid in non-stationary and non-vacuum spacetimes. We will argue that an appropriate 
definition of gravitational entropy should only involve the free gravitational field, as 
specified by the Weyl part of the curvature tensor, C a b cc i and that a particular 
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promising candidate for gravitational entropy can be constructed from integrals of 
quantities constructed from the pure Weyl form of the Bel- Robinson tensor |9] [10]. As we 
shall see below, this definition has the desired property of increasing as inhomogeneities 
form through gravitational attraction. It also reduces to the Bekenstein-Hawking value 
when evaluated in the case of a Schwarzschild black hole. 

Situations of particular interest are (i) those with Coulomb-like gravitational fields, 
representing the relativistic extension of Newtonian gravitational fields; and (ii) those 
involving wave-like gravitational fields. Each of these can be naturally represented by 
considering the Petrov classification of the Weyl tensor, and each constitutes a very 
different type of gravitational interaction. There is no a priori reason why one should 
expect the gravitational entropy in these two different settings to be describable in the 
same way, and we therefore consider them separately in what follows. As we shall see, 
this separation simplifies the problem of understanding the thermodynamic properties 
of each of these two types of gravitational fields considerably. The more general case, in 
which Coulomb and wave-like parts of the gravitational field are mixed, will, however, 
require a further and more detailed treatment. 

Throughout this paper we will make use of the 1 + 3 covariant description of 
gravitational fields [U [5], which proceeds by taking a timelike unit vector, u a , and 
defining a projection tensor, h a b = g a b + u a Ub- The covariant derivative of u a can then 
be split into irreducible parts such that 

U a -b = ~U a Ub + -Qh ab + (Tab + ^ab, (1) 

where = h ab u a -b is the expansion scalar, a ab = (h^h b d — \h a bh cd )u c ^d is the shear tensor, 
Wat = hr°h b fu c .d is the vorticity tensor, and u a = u b u a . b is the acceleration vector. The 
energy-momentum tensor of a fluid can then be decomposed relative to u a such that 

T ab = pu a u b + q a u b + u a q b + ph ab + rr ab , (2) 

where p is the energy density, q a is the momentum density, p is the isotropic pressure, and 
Ti a b is the anisotropic (tracefree) pressure of the fluid. In analogy to the decomposition 
of the Maxwell tensor into electric and magnetic parts, the Weyl tensor can also be 
decomposed into electric and magnetic parts as 

E ab = C ab cdU c u d and H ab = ^VacdC cd be u e , (3) 



where r] a bc = VabcdU d is the spatial alternating tensor (r] abcd = n [abcd]) 770123 = \f\g a b\), and 
both E ab and H ab are symmetric, tracefree and orthogonal to u a . We will also make use 
of the complex null tetrad 

m a = -i= (x a - iy a ) , l a = -i (u a - z a ) , and k a = 4= {u a + z a ) , (4) 

where x a , y a and z a are the spacelike unit vectors that, together with u a , form an 
orthonormal basis, and g a b = 2m^ a rhb) — 2k( a h) ■ 

In Section [2] we discuss our requirements for a sensible definition of gravitational 
entropy. These include being non-negative, vanishing when the Weyl tensor vanishes, 
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being a function of the anisotropy in the gravitational field, reproducing existing results 
for black hole entropy, and increasing monotonically in cosmological solutions when 
structure forms. In Section [3] we introduce and discuss the properties of the Bel- 
Robinson tensor, which plays a central role in our proposal. Section H] constructs a 
definition of entropy from the Bel-Robinson tensor, and from semi-classical notions 
of the temperature of a gravitational field. In Section [5] we apply our definition of 
entropy to the case of a Schwarzschild black hole, and find that it reproduces the 
Bekenstein-Hawking value. Section [6] then contains an application of our measure to 
scalar perturbations about a Robertson- Walker geometry, while in Section [7] we use the 
example of exact Lemaitre-Tolman-Bondi solutions to demonstrate the applicability 
of our proposed measure to non-perturbative inhomogeneous cosmological settings. 
Finally, in Section EJ we discuss our results. 

2. Requirements for Gravitational Entropy 

In general relativity the gravitational field can be split into Ricci and Weyl parts. The 
former is related pointwise to the energy-momentum tensor, where standard definitions 
for the entropy of matter fields should hold. Counting the entropy in the Ricci curvature 
of space-time would therefore be like counting the entropy in the matter fields twice. 
As our aim is to characterise the gravitational entropy of free gravitational fields, we 
therefore choose to concentrate on the Weyl part of the curvature tensor. This provides 
us with a tensorial description of the free part of the gravitational field, and is present 
even in the absense of matter fields. 

Definitions of gravitational entropy constructed from the Weyl tensor have 
been considered before, including the study of the simple choice S = C ab cd C cd ab . 
Unfortunately, this form of S has been argued to fail at isotropic singularities [6], and 
also fails to handle the decaying and growing perturbation modes [7]. Another choice 
has been the dimensionless scalar S = C ab cd C cd ab /R a b R b a , which while addressing the 
previous objections does not seem to give the correct sense of time for a radiating 
source [8], and also diverges for the vacuum case. Another shortcoming of these earlier 
definitions is that they are not always non-negative, and therefore cannot guarantee the 
monotonicity that is required of an entropy measure. They also appear as somewhat 
ad hoc, in that while they are motivated by the appropriate inclusion of the Weyl 
tensor in the numerator, they do not make any contact with either thermodynamics 
or information theory. The latter definition is also dimensionless, meaning it cannot 
reproduce the Bekenstein-Hawking result when applied to a black hole. 

To make progress on this issue we therefore make the following list of requirements 
on the gravitational entropy, S grav , that we expect to guide us: 

• El : It should be non-negative. 

• E2: It should vanish if and only if C a b c d = 0. 

• E3: It should measure the local anisotropy of the free gravitational field. 
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• E4 '■ It should reproduce the Bekenstein-Hawking entropy of a black hole. 

• E5: It should increase monotonically, as structure forms in the universe. 

As is the case for energy densities and pressures, we expect a sensible definition of 
gravitational entropy to be observer dependent (although it should be able to be 
defined covariantly) . While we are not considering quantum gravity here, it will 
also be considered beneficial if our definition of gravitational entropy can be linked 
to semi-classical calculations, in a similar way to the link between the Bekenstein- 
Hawking entropy and Hawking radiation. Finally, we note that we expect the entropy 
in gravitational and matter fields to be additive, such that the total of entropy in all 
fields is an extrinsic quantity. 

3. The Bel-Robinson Tensor 

With the above requirements in mind, we base our measure on the Bel-Robinson tensor, 
which is defined in terms of the Weyl tensor as [HI [TO] 

Tabcd = - (c eabf C e cd f + C* eabf C* e J j , (5) 

where C* abcd = ^r] abe fC e} \ d is the dual of the Weyl tensor. This tensor is overall 
symmetric, tracefree, and is covariantly conserved in vacuum (or in the presence of A). 
The factor of 1/4 here is included to give a natural interpretation of the Bel- Robinson 
tensor in terms of the Weyl spinor [TTJ. A measure of gravitational entropy constructed 
from this tensor has already been considered by Pelavas and Lake [12] and Pelavas and 
Coley [13] in the form S = J Wdr, where 

W = T abcd u a u b u c u d . (6) 

One can note immediately that W has the properties of being observer dependent and 
non-negative. We argue that observer dependence (i.e. dependence on u a ) is to be 
expected: this is the case for the energy density p and pressure p [3 0], and so may 
be expected also in this case. It also vanishes if and only if the Weyl tensor vanishes. 
This immediately addresses points El and E2. Similarly point E3 can be seen to be 
satisfied for perfect fluids, as the constraint equation for E ab in this case can be written 

D b E ab = -3co b H ab + ^D a p + [a, H] a , (7) 

where D a = h b Vh is the orthogonally projected 3-dimensional covariant derivative and 
[X, Y] a — rj abc X b d Y cd is a commutator [14] . Inhomogeneity in p then tells us that either 
E ab or H ab must be non-zero, and so we must have W > 0. Hence, inhomogeneity 
requires both anisotropy and a non-zero W. This is supported by the result that a 
conformally flat barotropic fluid solution is a Robertson- Walker geometry [15], which 
implies that W = if and only if the spacetime is homogeneous and isotropic at every 
point. 

Rather than constructing our entropy measure as an integral along a timelike curve, 
we shall employ integrals over spacelike hypersurfaces (see [IE]). We also wish to make 
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use of the gravito-electromagnetic properties of the Weyl tensor, as explored by Maartens 
and Bassett [14] . This starts with the recognition that the Bel- Robinson tensor given 
in Eq. (JH]) is the unique Maxwellian tensor [17] that can be constructed from the Weyl 
tensor, and that it acts as the "super-energy-momentum" tensor for gravitational fields 
[18], [19] . This can be seen by considering the 1+3 decomposition of the relevant evolution 
and constraint equations, which in vacuum are given by 0, H3] 

W=\ {E a b E\ + H a b H\) (8) 

J a = ~[E,H} a (9) 

W + D a J a ~ 0, (10) 

where J a = —h a e T e b C dU b u c u d is the "super-Poynting vector" and W is the "super-energy 
density", which can both be seen to arise naturally as invariants under spatial duality 
rotations in direct analogy with the energy density and Poynting vector of Maxwell's 
theory [14] . The ~ here indicates that terms that are expected to be small in the weak 
field limit have been discarded (for the full relativistic equation see [14]). The Bel- 
Robinson tensor therefore behaves like an effective energy-momentum tensor for the 
free gravitational fields, but has dimensions of L~ 4 , rather than the usual L~ 2 (where L 
is the unit of length). To find effective energy densities and pressures with the correct 
dimensionality we must therefore take a square-root. This problem has been considered 
by Bonilla and Senovilla [20], as we will now explain. 

For a symmetric and tracefree four-index tensor, such as T a b c d, one can define a 
symmetric two-index "square-root", t a b, as a solution of the following equation [22] : 

1 1 1 / 1 \ 

Tabcd = t{abtcd) — ^e(atb e 9cd) — ^t e et {ab9cd) + ^ ( t e ft B + j 9{ab9cd)-(^) 

The right-hand side of this equation is the only totally symmetric and tracefree four- 
index tensor that be constructed that is quadratic t a b- While it can be verified that 
for a given t a b this equation gives the only possible symmetric and tracefree four-index 
tensor, it is not the case that there exists a t a b for any arbitrary symmetric and tracefree 
four-index tensor. Furthermore, it can be seen that for any solution, t a b, there also 
exists another solution et a b + f g a b, where e = ±1 and where / is an arbitrary function 
[20] . We will explain in what follows how these ambiguities can be removed from our 
considerations for the cases that interest us here. 

The solution to Eq. ( 11 ip exists and is unique for tracefree t a b (i.e. when / = 0) in 
spacetimes of Petrov type D or N [20]. The tracefree requirement does not, however, 
necessarily lead to a quantity that is conserved in vacuum. One can therefore choose for 
the square-root of the Bel-Robinson to inherit its tracefree property, or (at least part 
of) its property of conservation in vacuum, but not necessarily both simultaneously. For 
the case of Petrov type D, in which there exist two double principal null directions, and 
the gravitational field is Coulomb-like, the tracefree square-root can be written as [20J 

t ab = 3e|^ 2 | (m( a rh b ) + l( a h)) , (12) 
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where ^2 — C a b c dk a m b rh c l d is the only non-zero Weyl scalar, and where we have used a 
complex null tetrad as defined in Eq. (j3J) with l a and k a aligned with the two principal 
null directions. For Petrov type N spacetimes, in which all four principal null directions 
are degenerate, and the gravitational field is wave-like, the tracefree square-root can be 
written as [20] 

tab = e\^ 4 \kah, (13) 

where = C a f )C drh a l b rh c l d is the only non-zero Weyl scalar in this case, and k a is chosen 
to be aligned with the principal null directions. This is obviously very similar to the 
energy-momentum tensor of pure radiation, and serves to further substantiate the claim 
that these objects behave like the effective energy-momentum of free gravitational fields. 

For other Petrov types the Bel- Robinson tensor can be factored into terms that are 
either Coulomb-like (as in the type D result above), wave- like (as in the case of type 
N), or more complicated. While the tracefree square- roots given above are unique, the 
factorisation of the Bel-Robinson tensor into two different symmetric tracefree tensors 
is not in general. This suggests that spacetimes for which the free gravitational fields 
are a mixture of wave-like and Coulomb-like parts could be more complicated, and that 
defining an effective energy density in free gravitational fields in these cases could also 
be complicated. Nevertheless, the cases given above contain very interesting examples 
for us to consider. These include all stationary black hole solutions, which are useful 
to compare to established definitions of gravitational thermodynamics, as well as the 
case of scalar perturbations to Robertson- Walker geometries, which are of great interest 
for cosmology. It is also the case that in all spacetimes with non-zero Weyl curvature 
a factor of the Bel-Robinson tensor exists which takes the form given in either Eq. 
(Fl2|) or ([13]) . This motivates using these expressions to determine the effective energy- 
momentum tensor for the Coulomb-like and wave-like parts of the free gravitational 
field, respectively. 



4. A Thermodynamic Description of the Free Gravitational Field 

In constructing our definition of gravitational entropy we will be seeking to construct a 
gravitational analogue of the fundamental law of thermodynamics in the form: 

-^grav^»Sgrav ^^grav Pgr&vdV, (14) 

where T grav , S gT!iv , ?7 grav and p gra v denote the effective temperature, entropy, internal 
energy and isotropic pressure of the free gravitational field, respectively, and V is the 
spatial volume. This will require a brief consideration of relativistic thermodynamics, as 
derived from the energy-momentum conservation equations, and the use of the square- 
root of the Bel- Robinson tensor to create an effective energy-momentum tensor for the 
free gravitational fieldjj]. 

\ Please note that by denning quantities such as C/ grav and p gra v we do not intend to imply that the 
free gravitational field should contribute to the right-hand side of the Einstein equations, only that 
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4-1. Relativistic Thermodynamics 

The subject of relativistic thermodynamics has been considered by a number of authors 
in the past (see, e.g., [2T] for a review). Starting from the equation for energy- momentum 
conservation we immediately find 

(u a T ab ) ]b = u a;b T ab -u a J a , (15) 

where J a = —T ab b , such that u a J a is heat flow into the fluid. If we now define = v/v, 
then Eqs. ©, © and (USD tell us that 

( pv y +pv = v [u a J a - q b . b - u a q a - a ab n ab ] , (16) 

where dots denote time derivatives along u a . This looks very much like a relativistic 
version of the fundamental thermodynamic equation, ( fl4l) . Given its form, it therefore 
seems natural to identify the left-hand side of this equation with Ts, where s and T are 
the point-wise entropy of the fields described by T ab , and their temperature, respectively. 
That is, we define 

Ts = (pv)' +pi), (17) 

where the entropy in a spatial domain a with volume V = J v will be given by S = f s. 
Clearly, in order to find S one first needs to define T independently. We shall return to 
this in Sec. 14.41 below. 



4-2. Effective Energy- Momentum of Coulomb-like Gravitational Fields 

Following the discussion in Section [3j we take the effective energy-momentum of the 
Coulomb-like gravitational fields that are present in a Petrov type D spacetime, Tab, to 
be given by the solution to Eq. ffTTj) . with a traceless part prescribed by Eq. f fl2|) . such 
that 

8nT ab = a [3e|* 2 | (m( a rh b) + l^k b) ) + fg ab ] (18) 
= a Qe|tf 2 | + /) (x a x b + y a y b ) - Qe|^ 2 | - f^j (z a z b - u a u b ) , 

where a is a constant to be determined, and in going to the second line we have 
transformed to a set of orthonormal basis vectors. We have used a curly T a b here 
to distinguish this quantity from the actual energy-momentum tensor of matter fields 
(as would appear on the right-hand side of Einstein's equations). 

We now contract this effective energy-momentum tensor with the unit 4-vector, 
u a , and the projection tensor, h ab , to gain the effective energy density, pressure and 
momentum density. These are given by 

871-^ =a 0e|* 2 | - A , 87rp grav = a0e|* 2 |+ A (19) 



47r7r o6 



rav a e\^ 2 \ (x a x b + y a y b - z a z b + u a u b ) , (20) 



quantities can be identified that appear to obey equations that are closely analogous to those of actual 
matter fields, and that these equations can be used to construct a definition of gravitational entropy 
that has many of the properties we require. 
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and gf rav = 0. The presence of the free function / in these expression is clearly 
undesirable, and can be removed by imposing the condition of energy conservation, 
u a T ab b = 0, in vacuum. As mentioned above, such a condition does not follow directly 
from the vacuum equation T abcd d = and Eq. (TTTj) . but must be imposed separately. 
In fact, we find it is satisfied if and only if 

/ = -i e |* 2 | + A 1 , (21) 

where Ax is an arbitrary constant that we will set to zero, as it does not affect the 
construction of the thermodynamic quantities that interest us. This result is found by 
differentiating Eq. ( Tig]) , and using |^ 2 | = y/2W/3 together with Eq. (45) from [Hj. 
The effective energy-momentum tensor is then given by 



2W 2W 
8nT ab = eaW — (x a x b + y a y b ) - 2ead — [z a z h - u a u ) , (22) 

and the effective energy density and pressure in the free gravitational field become 



WV av = 2aW— aIld Pgrav = 0, (23) 



'2W 

y 

where we have set e = +1, so that p grav > 0. The anisotropic pressure and momentum 
density are unchanged from Eq. (120!) and gf rav = 0, as / only occurs in the trace of 
Tab- These are the unique expressions for effective energy density and pressure of free 
gravitational fields in Petrov type D spacetimes that can be obtained from the square- 
root of the Bel- Robinson tensor under the condition of energy conservation in vacuum 
and positive energy densit}{§]. 

We can now construct from T a b in Eq. (|22|) an effective fundamental thermodynamic 
equation, of the form given in Eq. (fl6|) . In the presence of perfect fluid matter fields 
this is 



-^gravSgrav (Pgrav^O ^0~ab 



ab 47r(p + p) ■ 
gmv 3a^2W73 



(24) 



The first term in brackets on the right-hand side of this equation can be seen to be 
the direct analogue of the relativistic discipation that is present in Eq. (|16p . while the 
second term is due to the "heat flow" into the free gravitational fields from the matter 
fields. This term vanishes when p = p = 0, or when p = —p, or a a b = or E ab =0. While 
behaving as the analogue of heat in this equation, one should again bear in mind that 
the energy-momentum of the matter fields is being conserved as usual. This effective 
heat flow does not therefore signify any exchange of actual energy, but can be used in 
order to define entropy in the free gravitational fields as we have just done. 

To proceed further we need to quantify what we mean by gravitational temperature, 
T grav , in these settings. This we shall do in Sec. 14.41 Finally, let us note that the result 
|tf 2 1 = ^21^/3 means that the interesting properties of the Bel-Robinson tensor are 
inherited by our effective energy-momentum tensor T ab . This includes the satisfaction 
of points El, E2 and E3, as discussed in Section |3j 

§ Up to the possible inclusion of Ai, which would act like an effective vacuum term in Eq. ([22]) . 
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4-3. Effective Energy- Momentum of Wave-like Gravitational Fields 

Let us now consider plane-fronted transverse gravitational waves (of the sort that 
LIGO, and other gravitational wave experiments, are trying to detect). The geometry 
associated with these waves belongs to Kundt's class [23], which are the class of all Petrov 
type N solutions with vanishing Newman-Penrose scalar pnp = —k a]b m a fh b . These types 
of waves are closely analogous to our understanding of electromagnetic waves, and fit in 
with the idea of gravitational wave fronts as the congruence of null curves they follow 
is irrotational. 

Taking the effective energy-momentum tensor of this type of gravitational field to 
be given by the solution to Eq. ( TTTj) . with a traceless part prescribed by Eq. ( TlBl . leads 
to the effective energy-momentum tensor for the free gravitational fields being of the 
form 

^T ab = I3[e\^ 4 \k a k b + fg ab ] (25) 

where is a constant to be determined (which may or may not be equal to a in Eq. 
(JT8l)). In order to have energy conservation in vacuum in this case, such that u a T ab b = 0, 
we find that it is sufficient to take f — X2 =constant. That is, in this particular case, 
the tracefree condition and the energy conservation condition can both be satisfied 
simultaneously. 

The effective energy density, pressure, and momentum density of this effective fluid 
are then given by 

87rp grav =0 (||tt 4 | - A 2 ) (26) 
e 

e' 



87rp grav = /? ( ||^ 4 | + A 2 ) (27) 



8™ST = - P (|l*4| (x a x b + y a y b - 2z a z b )^ (28) 
8vrgf av = /^l^k, (29) 

and in the expressions for vr^ av and gf rav we have used a set of orthonormal basis 
vectors. We can now make the identification |\l/ 4 | = \J AW , and can consistently and 
without loss of generality set A2 = (as it is only the derivative of / that appears in 
the fundamental thermodynamic equation ffTB"]) ). Our thermodynamic quantities for the 
wave-like gravitational fields being considered in this section then become 

1 



87rp grav = (5\fAW and p grav = -p grav , (30) 

together with 7rf 6 av = -fiVAW (x a x b + y a y b - 2z a z b ) /48vr and 87rgf av = pVA~Wz a , 
where we have taken e = +1 so that p grav > 0. The effective fluid in this case therefore 
takes the form of radiation-like matter fields, with an equations of state w — p/p — 1/3. 
This can be compared with the dust-like equation of state that occurred in the case of 
Coulomb-like gravitational fields considered above. 
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Once again, from these considerations we can construct an effective gravitational 
analogue of the thermodynamic equation, f|T6|) . In the presence of a perfect fluid this is 

-^grav^grav (Pgrav^O Pgrnv^ 



pab 



= - v [<A-r + u a q r + a-] - 2nv(p ^ abEa • (si) 

Once more, the first term in brackets on the last line of this equation can be seen to be 
the direct analogue of the relativistic dissipation that is present in Eq. (1TB]) . while the 
second term is due to the "heat flow" into the free gravitational fields from the matter 
fields. This term again vanishes when p = p = 0, or p = —p, or when o a b = or E ab 
=0. For the case of general Petrov type N spacetimes there is an additional term on the 
right-hand side of this equation of the form 

4:71 



+ —VlWvk a ., b (x a x b + y a y h ). (32) 

While this term vanishes for the plane wave geometries in Kundt's class of solutions, 
it can be non-zero in general [21]. Finally, we note that the result |^4| = v4W again 
means that points El, E2 and E3 are satisfied. 

As an example of the energy-momentum of a type N solution one could consider 
an exact plane-fronted gravitational wave with line-element [25j [26] 

ds 2 = -dwdv + L 2 [e 2l dx 2 + e~ 2l dy 2 ] , (33) 

where L = L(w) and 7 = j(w). The vacuum field equations then reduce to 
L" + (7') 2 -^ = 0, where primes here denote differentiation with respect to w. The 
effective energy-momentum tensor given in Eq. (1251) can then be written as 

8nT ab = ^{iL 2 )'k a k b , (34) 

where k a = \/25 a v is the tangent vector to the single principal null direction. It is 
tempting at this point to compare this expression to the effective energy-momentum 
tensor for high-frequency gravitational waves found by Isaacson [27[ [28] . In this appoach 
one considers fluctuations of amplitude e < 1, and imposes the condition that derivatives 
acting on wave-like fluctuations are of order e _1 . For the solution above this corresponds 
to taking 7 ~ 0(e), so that L =constant+0(e). Eq. (131]) then becomes 

%TiT ab = 2(3 1 "k a k b + 0{l). (35) 

This expression is a factor of e _1 larger than the effective energy-momentum tensor 
found by Isaacson, which in the present case would be given by 

87r7 -Isaacson = 3 ( y } 2 ^ + (36) 

The effective energy-momentum tensor for wave-like gravitational fields that we have 
constructed above is therefore not the same as the one that appears on the right-hand 
side of the field equations in Isaacson's approach. This is not too surprising, as we have 
constructed our definition from the Weyl curvature, while Isaacson's definition is taken 
to be proportional to the Ricci curvature through the field equations. 
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4-4- Temperature of Gravitational Fields 

In the discussion presented so far in this section we have calculated expressions for 
TgraySgj-av by appealing to the analogy between the effective energy-momentum tensor 
of the gravitational fields that we have constructed, and the actual energy-momentum 
tensor of matter fields. In order to go from these equations to a calculation of entropy 
specifically we now need to know the "temperature", T grav , of the free gravitational 
fields. It is normally the case that, when considering the thermodynamics of ordinary 
matter fields, it is not possible to determine the temperature of a fluid without knowing 
something about the underlying microscopic physics. In the kinetic theory of gases, 
for example, one can identify temperature with the average kinetic energy of the gas 
molecules, but one cannot determine temperature from consideration of the macroscopic 
thermodynamic variables only. In this regard, it is not unreasonable to assume that a 
thermodynamic consideration of the free gravitational fields should be any different to 
that of standard thermodynamics: In order to know the temperature one needs to know 
something of the physics of the underlying microscopic theory. 

For the temperature of gravitational fields we therefore appeal to the results of 
black hole thermodynamics pQ, and quantum field theory in curved spacetimes [29]. For 
our purposes we require a definition of temperature that is local (rather than being 
defined for horizons only), and that reproduces the expected results from semi-classical 
calculations in spacetimes such as Schwarzschild and de Sitter. We therefore take the 
temperature at any point in spacetime to be given by the following expression: 

\u a , b l a k h \ _ \u a z a + H + a ab z a z b \ 

Jgrav ^ 2n , 

where l a = (u a — z a )/V2 and k a = (u a + z a )/V2 are the real vectors in a complex null 
tetrad, z a is a spacelike unit vector aligned with the Weyl principal tetrad, and H = 6/3 
is the isotropic Hubble rate. This expression reproduces the Hawking temperature 
[TJ, the Unruh temperature [30] and the temperature of de Sitter space pH] in the 
appropriate limits. One should note that Eq. fl37|) is an extra ingredient in our analysis, 
beyond the construction of Tab, an d one could use alternative definitions as appropriate^. 



5. The Entropy of Black Holes 

The entropy of stationary black hole spacetimes is of obvious importance in discussing 
gravitational entropy, as they are so far the only spacetimes that allow an unambiguous 
definition of entropy in free gravitational fields. We will therefore consider them here, in 
the light of our previous discussion, in order to investigate whether our proposal agrees 
with previously established results. 

|| This may, for example, be the case if one wishes to try and use perturbed Friedmann-Lemaitre- 
Robertson- Walker solutions as effective macroscopic descriptions of more complicated microscopic 
spacetimes. The underlying theory would then be Einstein's theory, and the temperature should 
presumably arise out of the consideration of averaging or coarse-graining procedures. 
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We consider the Schwarzschild geometry written in Gullstrand-Painleve coordinates 

as 



als 2 = -(l-^j dt 2 - 2^^drdt + air 2 + rW, (38) 

where m is the constant mass parameter. This coordinate system admits hypersurfaces 
of constant t that intersect the horizon, and that have Euclidean geometry. For our unit 
vectors u a and z a we take 

« a = [ 0;^^=,0,0 | (39) 

j 2m I 



0,0,0 , (40) 



2m I 



such that in the region r < 2m we have u a u a = — 1 and z a z a = 1. It is also the case 
that u a z a = 0, and z a is orthogonal to the Euclidean hypersurfaces with t =constant. 
This choice is such that u a and z a specify a Weyl principal tetrad, as in Eq. (@J. 

From Eqs. ( 123]) and ( 137|) we then have that the gravitational energy density and 
temperature is given at each point in the region r < 2m by 

2m 

8vrp grav = a— (41) 

Tgrav = / =• (42) 

27rrV|l - 22* I 



Now, Eq. f )24p cannot be directly applied to stationary spacetimes as strictly it involves 
changes in thermodynamic quantities over time. Following the methods used in [TJ, we 
therefore instead choose to compare two different stationary black hole spacetimes that 
might represent the gravitational field before and after a small amount of mass is added 
to the black hole. Eq. ( 1241) then becomes 

iV = %=^, (43) 

^ grav 

for an incremental increase in effective gravitational energy at a given T grav along each 
of the curves u a . Integrating over a volume V on a hypersurface of constant t then gives 

5 grav = / Vav= / [ ^wy6drd0d4, (44) 

JV JV -'grav Jv ^ 

where we have taken v = z a 7] a b c ddx b dx c dx d and set an arbitrary constant to zero. If V 
is the region interior to the event horizon then this gives 

Sgrav = tt^p, (45) 

where = 4ir(2m) 2 is the area of the horizon. This expression is identical to the 
Bekenstein-Hawking value if a — 1, and therefore satisfies point E4. 

In order to consider the entropy in more general black hole spacetimes it will 
be necessary to include the entropy associated with the extra degrees of freedom 
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involved. For the Reissner-Nordstrom solution this will mean including the energy 
density associated with the electromagnetic field, while in the Kerr solution it will 
require taking into account the energy associated with the rotation. Nevertheless, we 
see no reason in principle why the above approach could not be applied to these solutions. 
They are all Petrov type D, so the same equations should be valid. We leave this for 
future work. 



6. The Entropy of Large-Scale Structure in Cosmology 

Let us now consider a spatially flat Robert son- Walker geometry with scalar 
perturbations in a longitudinal gauge, such that the line-element can be written 

ds 2 = a 2 (r) [-(1 + 2<p)dT 2 + (1 - 2<p)(dx 2 + dy 2 + dz 2 )] . (46) 

Our reference set of curves in this spacetime will be given by 

ti a = (47) 
1 (O;V,0), (48) 



a|V0| 

where i,j are spatial indices, and higher-order terms have been discarded. 

To lowest order in perturbations we then have that the space-time is silent, with 
H a b = 0. Meanwhile, the electric part of the Weyl tensor is given in tetrad components 
by 

3 

E AB = e A a e B b E ab = (f). AB + J2 ( t>.AlCAB for A^B (49) 

c=i 

where ■jabc = e J 4 a; f ) e B a e c 6 are the Ricci rotation coefficients, and a dot denotes the 
tetrad component of a covariant derivative such that 4>.ab = e B 6 (e yl a ;a ) ; b. By aligning 
e3 with V0 we have <p,i = = <p,2 everywhere, and it can be deduced from Eq. ( 1491) that 
in this case E ab is diagonalized with E ab = Y1a=i ^A^Aa^Ab, where A^ are the eigenvalues 
of this tensor. The Petrov type of the spacetime in Eq. ( 1461) must therefore be either 
I, D or O, to lowest order in perturbation^. That the diagonal form of E ab should be 
maintained after spatial rotations means that Ai = A2, so the spacetime is in fact type 
D (unless (ft =constant, in which case it is type O). 

Aligning z a with e 3 , as in Eq. (JIB"]) , and using Eqs. (123|) . ( 1371) and (j47j) . we find the 
gravitational energy density and temperature to be given to lowest order by 

. . \{a A uu z l z : >\ 

8vrp grav = a I <p,{ij)Z z J \ = a- — (50) 

a 

^grav = ^ZT' (51) 

where a dot denotes a derivative with respect to r, and angled brackets around indices 
denote that a quantity is tracefree. 



If What we mean by this is that the lowest order scalar perturbations in Eq. (|46|) lead to a Weyl tensor 
that is of one of these Petrov types, when all higher order terms are neglected. 
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After integrating over a spatial volume V, Eq. f )24p now gives the evolution of the 
gravitational entropy of this spacetime to lowest order as 

4- = ^H\j v ( ai y/Wfi) dxdydz (52) 

j dxdydz, (53) 

where we have used v = a?dxdydz. We take the spatial domain V to be a comoving 
spatial volume, although one is free to make other choices. It can be seen that our choice 
of temperature in Eq. (jSTl) reproduces the expected Gibbons- Hawking temperature of 
de Sitter space [31J, as well as the temperature of the horizon in other homogeneous and 
isotropic spaces [32J. 

This formulation of gravitational entropy can be seen to be directly dependent 
on the shear of u l , and on the anisotropy of the gravitational field through It 
therefore explicitly satisfies point E3. Factoring and u l into time-dependent and time- 
independent parts, as is usual in the study of perturbed Robertson- Walker geometries, 
it can be seen that the approximate time dependence of S grs , v goes like 

Sgrav ~ a3 U, (54) 

where ui = u(r)Qi(x^), and where Qi{x^) is a harmonic function. Now for the 
monotonicity condition E5 to be satisfied, the quantity on the right-hand side of Eq. fl54|) 
needs to be a monotonically increasing function of time. For expanding dust dominated 
universes we note that Eq. ( 1541 gives 

S grav ~ r 5 ~ t 5 / 3 , (55) 

where t — J a(r)dr is the proper time of comoving observers. This certainly satisfies the 
monotonicity condition. In fact, the S gTav in Eq. (l54"j) will grow for any u that decays 
more slowly than ~ a -3 . 



a 
4W\ 



\ 



(a u 



) z l z 3 



7. A Non-perturbative Inhomogeneous Example 

It would be of interest to see whether the above measure of entropy is also applicable 
in non-perturbative inhomogeneous settings. To this end, let us consider the Lemaitre- 
Tolman-Bondi (LTB) solution, with line-element given by 

ds 2 = -dt 2 + - % r + R 2 dQ 2 , (56) 
1 — k[r) 

where R = R(r, t), and prime denotes partial differentiation with respect to r. The field 
equations then reduce to 

R 2 m(r) k(r) m! 

lP = ^r-lv and p = ^Rm>' (57) 

where dots here denote differentiation with resect to t. The free function m(r) in this 
solution specifies a measure of the gravitational mass within a sphere of radius r, and 
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in the FLRW limit reduces to m — > 4ttR 3 p/3. The LTB spacetime is of Petrov type D, 
unless it is FLRW. 

If we now take our time-like and space-like unit vectors to be 

u a = (1; 0,0,0) (58) 

*•= (o; S,o, oV (59) 

then we obtain the following expressions for the effective energy density and temperature 
of the free gravitational fields 

87rp grav = 2a l - R3 H/ 1 (60) 



T = ±- 

ZTT 



R' 



R' 



(61) 



The effective energy density grows as the departure from FLRW increases, weighted by 
a factor of R 3 , and the effective temperature is simply proportional to the expansion in 
the radial direction. 

In a given LTB spacetime, entropy will then increase as long as the expansion 
weighted value of p gra v increases. For a given LTB solution, we therefore have an arrow 
of time that associates low p gra v with early times, and high p grav with late times. This 
provides us with a way of specifying which direction is "the future" for any given solution 
(independent of whether the coordinate t happens to be increasing in that direction or 
not, or whether the spacetime is expanding or not) [33]. It also supports our intuitive 
understanding of gravitational entropy in cosmology, in which the low entropy state 
of the early universe was close to homogeneous and isotropic, and in which the late 
universe is inhomogeneous and anisotropic. 

What this means for an example such as a time-reversal-symmetric recollapsing 
LTB solution that starts off close to homogeneous, then evolves towards inhomogeneity 
at its maximum of expansion, before recollapsing to a near homogeneous final state, 
remains to be seen. It seems safe to conjecture that solutions of this type will always be 
problematic for any sensible definition of gravitational entropy, however, as time reversal 
symmetry about the maximum of expansion will necessarily mean that any notion of 
entropy based on the geometry of spacetime will not be monotonic throughout its entire 
evolution. An example such as this would appear to have an entropic arrow of time 
that points towards the maximum of expansion from both sides, so that the "future" 
(defined in this way) is always in the direction of positive cosmological expansion. 



8. Discussion 



Motivated by thermodynamical considerations, we have proposed a measure of 
gravitational entropy based on the square-root of the Bel-Robinson tensor. A key 
feature of this measure is its non-negativeness, which is a fundamental requirement 
for any measure of entropy, and is in contrast to other gravitational entropy measures 
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previously proposed. We have applied our measure to a number of examples, including 
cosmological ones, and have found that the specific form of this measure will depend on 
whether the gravitational field is Coulomb-like or wave-like. 

In constructing this measure of gravitational entropy we have assumed that the 
square-root of the Bel-Robinson tensor can be taken to be the effective energy- 
momentum tensor of free gravitational fields, that the energy of this effective 
fluid is conserved in vacuum, that the temperature of gravitational fields can be 
determined from semi-classical considerations, and that the fundamental equation of 
thermodynamics is applicable. We have also assumed that the entropy of the universe at 
any given time is given by integrating the entropy density over a space-like hypersurface. 
Under these assumptions our definition of entropy in free gravitational fields is unique. 
For the square-root to exist as a unique factorisation of the Bel-Robinson, however, 
requires the spacetime to contain gravitational fields that are only Coulomb-like or only 
wave-like. The factorisations of the Bel-Robinson tensor that are possible when both 
of these types of field are present are more complicated, and are not unique. Further 
study is therefore necessary to extend this definition to the case of general spacetimes. 

In the case of the Schwarzschild black hole, our entropy measure reduces to the usual 
Bekenstein-Hawking value. For scalar perturbations of a Robertson- Walker geometry 
we find that our measure evolves like the Hubble weighted anisotropy of the free 
gravitational field. As a result it increases as structure formation occurs, as is expected 
of a sensible measure. We have also applied our measure to the non-perturbative case of 
LTB models and found conditions under which the entropy increases. These examples, 
ranging from black holes to perturbed FLRW and exact inhomogeneous LTB, provide 
encouraging evidence that our proposed measure has the potential of accounting for 
gravitational entropy in a range of cosmological settings of interest. However, it is only 
in the black hole case that there is and obvious link to the holographic principle [33] . 

Finally, we note that the definition of gravitational entropy we have considered in 
this paper is only valid for General Relativity. In other theories of gravity the Bel- 
Robinson may not be the appropriate choice to describe the super-energy-momentum 
of free gravitational fields, and a Bel- Robinson-like tensor may not exist at all [35J. 
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